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MOLECULAR PHYLOGENY OF THE EDELWEISS 
(LEONTOPODIUM, ASTERACEAE - 
GNAPHALIEAE) 


C. BrócH!, W. B. Dickon£É?, R. SAMUEL! & T. F. STUESsyY! 


Leontopodium is a genus of approximately 30 species with a conspicuous Asian-European 
disjunct distribution. In this study samples from the Himalayan/Tibetan centre of 
diversity of the genus, as well as from Europe, were analysed to infer a phylogeny of the 
genus using sequences of nuclear ribosomal (ITS and ETS) and plastid (matK and trnL/F) 
DNA. The Southeast Tibetan monotypic Sinoleontopodium [Leontopodium lingianum 
(Y.L.Chen) Dickoré, comb. nov.] falls into Leontopodium. Monophyly of Leontopodium, 
including Sinoleontopodium, is supported. Due to low rates of sequence divergence, 
intrageneric relationships in general are weakly supported, a pattern frequently observed 
in plant groups centred in the Tibetan Plateau. Three phylogenetic groups can be 
identified, however, and these are also supported by morphology. The low levels of 
nucleotide divergence suggest a young age for the group, which has been influenced by the 
turbulent geological history of the Tibetan Plateau. Leontopodium is a characteristic Sino- 
Himalayan element that appears to have found its way into the mountains of Europe in 
geologically recent times. The two European taxa, Leontopodium alpinum and L. nivale, 
form a genetically distinct group, which, considering the wide geographic disjunction, 
shows surprisingly little divergence from its Asian relatives. 


Keywords. Asteraceae, ETS, Gnaphalieae, ITS, Leontopodium, matK, molecular phyl- 
ogeny, Sino-Himalayan flora, Sinoleontopodium, trnL/F. 


INTRODUCTION 


The alpine edelweiss (Leontopodium alpinum Cass. or L. nivale subsp. alpinum (Cass.) 
Greuter*) holds an important position in the cultural heritage of the European Alps. 
Starting from a tradition in folk medicine, phytochemical research has revealed the 
presence of various new or uncommon secondary compounds, some of which are of 
high biological activity (Grey et al., 1999; Stuppner et al., 2002; Dobner et al., 2003a, 
2003b, 2004; Schwaiger et al., 2004). The present study has evolved under the growing 
need for an applicable taxonomy and systematics of the genus for further pharmaceu- 
tical research. 
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Leontopodium (Pers.) R.Br. is a genus containing about 120 published binomials. 
Based on the earlier work of Franchet (1892) and Beauverd (1909, 1910, 1911, 1912, 
1914), Handel-Mazzetti (1927) completed the only monograph to date, in which he 
accepted 41 species. Several new taxa have subsequently been described (e.g. 
Voroshilov, 1979; Akiyama, 1999). We found the species concepts of Handel-Mazzetti 
(1927) generally sound. However, several taxa he accepted, and almost all of those 
described since, are probably better reduced to synonymy. Many were known only 
from their types and, on account of more recent collections, are now known to fall 
within the often wide morphological variation of other species. Ongoing studies 
suggest there are fewer species, perhaps as few as 30 (W. B. Dickore, unpubl.). 

Most taxonomic treatments of Leontopodium in modern Asian Floras, such as the 
Flora of the USSR (Shishkin, 1999), Flora Iranica (Rechinger, 1980), Flora of India 
(Hajra et al., 1995), Flora of Pakistan (Qaiser & Abid, 2003), Flora of Japan 
(Iwatsuki et al., 1995), Flora of Bhutan (Grierson & Long, 2001) and for Central Asia 
(Grubov, 2003), do not seem fully satisfactory, because of their rather incompatible 
regional viewpoints. 

Except for two taxa in the mountains of Central and Southern Europe 
(Leontopodium alpinum and L. nivale; author abbreviations for species are given in 
Table 1), the genus occupies an extensive and largely continuous distribution in the 
mountains of Central, temperate Southeastern, and Eastern Asia. This wide generic 
disjunction is unusual (Meusel & Jager, 1992). Whereas Leontopodium alpinum is 
well known due to ample work on various aspects of its biology, ecology, chorology, 
and secondary chemistry (e.g. Uexküll-Gyllenband, 1901; Sokolowska-Kulczycka, 
1959; Maugini, 1962; Tira et al., 1970; Erhardt, 1993; Hook, 1994; Stuppner et al., 
2002), very little is known about the Asian members of the genus. 

This study employs DNA sequences to understand affinities among species of 
Leontopodium and among some possible close generic relatives. More than two- 
thirds of the species of Leontopodium are sampled, hence providing a good 
taxonomic and geographical representation, plus taxa from related genera of 
Gnaphalieae and more distantly related taxa from other tribes of Asteraceae. Two 
nuclear ribosomal DNA markers (internal transcribed spacer (ITS) and external 
transcribed spacer (ETS)) and three plastid markers (matK, trnL intron, and trnL/F 
intergenic spacer) were investigated. ITS1 and ITS2 have frequently been used for 
phylogenetic studies in Asteraceae especially at the generic level (Clevinger & 
Panero, 2000; Samuel et al., 2003, 2006; Kimball & Crawford, 2004; Plovanich & 
Panero, 2004; Tremetsberger et al., 2005), and for phylogenetic studies in Gnapha- 
lieae (Antennaria Gaertn., Bayer et al., 1996; Anaphalioides (Benth.) Kirp., Glenny, 
1997; Australasian Gnaphalieae, Breitwieser et al., 1999; Mediterranean Helichrysum 
Mill, Galbany-Casals et al., 2004a). The ETS region of the 35S rDNA cistron 
belongs to the same functional unit as ITS (Baldwin & Markos, 1998), but so far 
ETS has been used less frequently (e.g. Baldwin & Markos, 1998; Clevinger & 
Panero, 2000; Plovanich & Panero, 2004). The matK gene, which is one of the fastest 


TABLE 1. Species names, localities, collectors and voucher numbers, and GenBank accession numbers of the analysed taxa. All vouchers are 
deposited in WU unless stated otherwise. Collectors are abbreviated as follows: AT = A. Tribsch; BD = B. Dickoré; EH = E. Hörandl; FE = 
F. Essl; GM = G. Miehe; GS = G. Schneeweiss; I = Iwatsuki; JR = J. F. Rock; K = Kleesadl; KR = K. Reiter; LK = L. Klimes; MD = M. 
Dobner; MW = M. Wiedermann; N = Narantuja; PS = P. Schönswetter; SE = S. Ertl; SM = S. Miehe; S = Smith; UW = U. Wündisch; 
WT = W. Till; Y = Yamazaki. Voucher information for sequences drawn from GenBank can be found in the indicated references 


Species 


Leontopodium 
alpinum Cass. 


Accession Locality 


1 


2 
3 


Italy, Piemonte 


Spain, Aragon 

Spain, Aragon 

Italy, 
Lombardia, 
Sondrio 


France, Alpes de 
Haute Provence 


Italy, Alto 
Aldige 
France, Savoie 


Italy, Lombardia 


Italy, Cuneo 

Austria, 
Kärnten 

Italy, Trento 


Collector, 
collection 
number and 
herbarium 
acronym 


GS & PS 8914 
GS & PS 8959 
GS & PS 8881 
MW 9282 
MW 9029 

PS & AT 9016 
SE 9377 

GS & PS 8940 
MW 9197 

PS 9059 


GS & PS 8809 


ITS 
FJ39910 
FJ639911 
FJ639912 
FJ639913 
FJ639914 
FJ639915 
FJ639916 
FJ639917 
FJ639918 
FJ639919 


FJ639920 


ETS 


FJ639981 


matK. 


FJ640024 


trnL 
intron 


trnL/F 


AF141821 


AF141733 — 


Reference 


This study 


This study 
This study 
This study 


This study 
This study 


This study 
This study 
This study 
This study 


This study 
Bayer 

et al., 2002 
Bayer 

et al., 2002 
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Species 


L. andersonii 
C.B.Clarke 


L. artemisiifolium 
(H.Lev.) Beauverd 


L. caespitosum Diels 


L. calocephalum 
(Franch.) 
Beauverd 


L. dedekensii 
(Bureau & 


Franch.) Beauverd 


Accession Locality 


1 


2 
3 


China, Yunnan 


China, Yunnan 
China, Yunnan 


China, Yunnan 
China, Yunnan 


China, Yunnan 
China, Yunnan 


China, Yunnan 
China, Yunnan 


China, Yunnan 
China, Sichuan 
China, Yunnan 


China, Yunnan 


China, 
Tibet/Xizang 


Collector, 
collection 
number and 
herbarium 
acronym 


MD et al. 
MD201-2394 
BD 14068 
BD 14570 
(B, WU) 
I et al. s.n. 
(BM) 
MD et al. 
D01-2396 
BD 14574 
MD et al. 
D01-2384 
BD 14040 
MD et al. 
MD201-2354 
MD et al. 
MD201-2395 





UW 94-488-5 
MD et al. 
MD201-2358 


BD 14159 
(B, WU) 
BD 8664 (B) 


ITS 
FJ639921 


FJ639922 


FJ639923 


FJ639924 
FJ639925 


FJ639926 
FJ639927 


FJ639929 


FJ639928 
FJ639930 


FJ639932 


FJ639931 


ETS 


FJ640006 
FJ640007 


FJ640008 


FJ640009 
FJ639982 


FJ640010 
FJ640011 


FJ640012 


FJ639994 
FJ640001 


FJ640002 


matK. 


FJ640038 


FJ640025 


FJ640026 


trnL 
intron 


FJ640041 


FJ640042 
FJ640043 


trnL/F 


FJ639971 


FJ639972 
FJ639973 


Reference 


This study 


This study 
This study 


This study 
This study 


This study 
This study 


This study 
This study 


This study 


This study 
This study 


This study 


This study 
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L. discolor 
Beauverd 


L. franchetii 
Beauverd 


L. haastioides 
Hand.-Mazz. 


L. himalayanum DC. 


L. jacotianum 
Beauverd 


L. japonicum Mig. 
L. leontopodinum 
(DC.) Hand.-Mazz. 


L. leontopodioides 
(Willd.) Beauverd 


China, Yunnan 


China, Yunnan 
China, Yunnan 
China, 
Tibet/Xizang 
China, 
Tibet/Xizang 
China, 
Tibet/Xizang 
China, 
Tibet/Xizang 


Nepal, Langtang 


Nepal, Langtang 


China, Sichuan 


Nepal, Langtang 


China, 
Tibet/Xizang 

Japan, Shinano 

India, Jammu 
& Kashmir 


China, Qinghai 


China, Gansu 
Mongolia, 


Hövgöl Aimak 


MD et al. 
MD01-2411 
BD 14499 
BD 14495 
BD 9538 (B) 


BD 9892 
(B, WU) 
BD 5228 
(B, WU) 
BD 9945 (B) 


KR 571 


KR 25 

MD et al. 
MD01-2445 

KR518 

BD 9167 (B) 


T 1228 (W) 

Southampton 
University 
51 (BM) 

GM & SM 
9315/12 (B) 

JR 12239 (BM) 

N S-070800 
(BM) 


AF002168 — 


FJ639933 — 
FJ639934 FJ639983 
FJ639935 FJ640003 
FJ639936 — 
FJ639937 FJ640013 
FJ639938 FJ640004 
FJ639939 — 


FJ63994] — 


FJ639940 FJ639995 


FJ639942 — 
= FJ64001 


= FJ639984 
= FJ639987 
FJ639943 — 


- FJ639986 
FJ639944 FJ640014 


FJ640027 


FJ640044 


FJ640046 


FJ640045 


FJ639974 


FJ639976 


FJ639975 


Glenny & 
Wagstaff, 
1997 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 


This study 
This study 


This study 
This study 


This study 
This study 
This study 


This study 
This study 
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Species 


L. microphyllum 
Hayata 


L. nanum 
Hand.-Mazz. 


L. nivale (Ten.) Huet 
ex Hand.-Mazz. 


L. ochroleucum 
Beauverd 


L. pusillum 
(Beauverd) 
Hand.-Mazz. 


Accession Locality 


1 


2 
3 


Taiwan, Hualien 


Taiwan, Hualien 
Taiwan, 
Nanhuta shan 
China, Xizang, 
Nagarze, 
Pomo Co. 
China, Qinghai 


Italy, Abruzzo 


Italy, Abruzzo 
India, Jammu & 
Kashmir 
Russia, Altay 
Russia, Altay 
China, 
Tibet/Xizang 


China, Qinghai 


China, Qinghai 


Collector, 

collection 

number and 
herbarium 

acronym ITS 


EH et al. 9549 | FJ639946 


EH et al. 9555 FJ639947 
Y et al. 303 (BM) — 


BD 9539 (B) FJ639951 

GM et al. FJ639948 
98-34117 (B) 

GS & PS 8926 FJ639950 


EH et al. 6244 | FJ639949 
LK 03-21-30 FJ639954 


AT 9523 FJ639952 
AT 9594 FJ639953 
GM et al. FJ639955 


97-073-06 (B) 


BD 4425 FJ639956 
(B, WU) 

GM et al. — 
98-35212 (B) 


ETS 
FJ640015 


FJ640016 
FJ639989 


FJ640005 


FJ640017 
FJ639997 


FJ639996 
FJ639990 


FJ640018 


mark 


FJ640040 


FJ640039 


FJ640029 


FJ640030 
FJ640032 


FJ640031 
FJ640033 


FJ640034 


trnL 
intron 


FJ640047 


FJ640048 


trnL/F 
FJ639977 


FJ639978 


Reference 


This study 


This study 
This study 


This study 


This study 
This study 


This study 
This study 


This study 
This study 
This study 


This study 


This study 
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L. sinense Hemsl. 


L. souliei Beauverd 


L. stracheyi 
(Hook.f.) 
C.B.Clarke 
ex Hemsl. 


L. cf. stracheyi 
(Hook.f.) 
C.B.Clarke ex 
Hemsl. 

Sinoleontopodium 
lingianum 
Y.L.Chen 


OUTGROUPS 
Achillea 
millefolium L. 


Anaphalioides 
subrigida 


(Colenso) Anderb. 


China, Yunnan 


China, Yunnan 
China, Sichuan 
China, Yunnan 


China, 
Tibet/Xizang 


China, 
Tibet/Xizang 

China, 
Tibet/Xizang 


China, 
Tibet/Xizang 


China, 
Tibet/Xizang 


MD et al. 
MD201-2397 
BD 14534 (B) 
S 12219 (BM) 
MD et al. 
MD01-2404 
GM & SM 
98-09509 (B) 


BD 3404 (B) 


BD 10529 (B) 


BD 10836 


BD 11363 
(B, WU) 


FJ639958 
FJ639957 
FJ639959 


FJ639960 


FJ639961 


FJ639945 


U95279 


FJ639999 — 
FJ640000 FJ640035 
FJ63999] — 
FJ639998 FJ640036 


FJ640020 FJ640037 


FJ640021 — 


FJ640022 FJ640028 


FJ639988 — 


= EU385315 - - 


This study 


- - This study 
This study 


FJ640049 FJ639979 This study 


This study 


This study 


This study 


This study 
This study 


EU385030 EU385030 Panero & 
Funk, 2008 

Panero & 
Funk, 2008 

Glenny & 
Wagstaff, 
1997 
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TABLE l. (Cont’d) 


Species Accession Locality 


Anaphalis javanica 
Sch.Bip. in Zoll. 


Anaphalis 
margaritacea 
Benth. & Hook.f. 

Anaphalis nubigena 
DC. 


Antennaria 
anaphaloides 
Rydb. 

Antennaria carpatica 
(Wahlenb.) Bluff & 
Fingerh. 

Antennaria dioica (L.) 1 
Gaertn. 

2 


Antennaria luzuloides 
Torr. & Gray 
Castroviejoa frigida 
(Labill.) Galbany, 
L.Sáez & 
Benedi 


India, Jammu & 


Kashmir, 
Ladakh 


Austria, 
Salzburg 


Austria, 


Niederósterreich 


Russia, 
Altay 


Collector, 
collection 
number and 
herbarium 
acronym 


LK 04-01-9 


EH H63-30.8.88 


WT s.n. 


AT & FE 
004588-8978 


ITS ETS 


U95292  - 


AF046937 — 


FJ639962 FJ640023 


L40857 - 


FJ639965 — 


FJ639964 FJ639985 


= FJ639992 


AY445228 — 


trnL 
matK intron 


AF456774 — 


FJ640050 


trnL/F 


FJ639980 


Reference 


Glenny & 
Wagstaff, 
1997 

Noyes & 
Rieseberg, 
1999 

This study 


Bayer et al., 
1996 


This study 


This study 
This study 


Bayer et al., 
2000 


Galbany-Casals 


et al., 2004a 
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Castroviejoa 
montelinasana 
(Em.Schmid) 
Galbany, 

L.Sáez & 
Benedi 

Ewartia argentifolia 
N.A.Wakef. 

Ewartia catipes 
(Hook.f.) Beauverd 


Ewartia meredithiae 
(F.Muell.) Beauverd 

Ewartia planchonii 
(Hook.f.) Beauverd 

Ewartiothamnus 
sinclairii (Hook.f.) 
Anderb. 

Filago pyramidata L. 


Gamochaeta 
pensylvanica 
(Willd.) 
Cabrera 


1 


2 


3 


AY445229 — 


AF115910 — 


= FJ404694 


U95290 — 


AF115907 — 


AF115909 — 


U95283 - 


AY445231 — 


= - = Galbany-Casals 
et al., 2004a 


— — — Breitwieser 
et al., 1999 

AF151460 AF141698 AF141786 Bayer et al., 

2002 
R. D. Smissen, 
unpubl. 

- -— — Glenny & 
Wagstaff, 
1997 

= = = Breitwieser 
et al., 1999 

= = = Breitwieser 
et al., 1999 

- - - Glenny & 
Wagstaff, 
1997 

- - - Galbany- 
Casals 
et al., 
2004a 

= EU385070 EU385070 Panero 
& Funk, 
2008 
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Species 


Gamochaeta spicata 
(Lam.) Cabrera 
Gnaphalium 
hoppeanum 
(L.) Hilliard & 
B.L.Burtt 
Gnaphalium 
norvegicum 
Gunnerus 
Gnaphalium 
supinum L. 


Gnaphalium 
sylvaticum L. 
Helichrysum 


leucopsideum DC. 


Leucogenes 
grandiceps 
(Hook.f.) 
Beauverd 


Accession Locality 


Italy, 


Lombardia 


Austria, 


Oberösterreich 


Austria, 


Niederösterreich 


1 New Zealand, 
Canterbury 


Collector, 
collection 
number and 
herbarium 
acronym 


PS & AT 5382 


WT s.n. 


PS et al. 
004137-8982a 


ITS ETS 
AF115917 — 


FJ639967 — 


FJ639968 — 


AY445230 — 


FJ639969 FJ639993 


matK. 


- AF319712 — 


FJ639970 — 


trnL 
intron 


trnL/F 


Reference 


Breitwieser 
et al., 1999 
This study 


This study 


Galbany- 
Casals 
et al., 
2004a 

This study 


Bayer 
et al., 2002 
This study 


AY606885 AY 606896 Smissen et al., 


2004 





ve 
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TABLE 1. (Cont'd) 


Leucogenes 
leontopodium 
(Hook.f.) Beauverd 

Leysera gnaphalodes 
(L.) L. 


Pseudognaphalium 
luteoalbum 
(L.) Hilliard & 
B.L.Burtt 
Oedera squarrosa 
(L) 
Anderb. & 
K.Bremer 


Relhania calycina 
L.Hér. 

Solidago rugosa 
Mill. 

Symphyotrichum 
novae-angliae 
(L.) G.L.Nesom 


AF543357 — 


_ D = Smissen et al., 
2003 


= — AF100473 Bayer et al., 
2000 

AF456793 — = Bayer & 
Cross, 2003 

- - - Galbany- 
Casals 
et al., 
2004a 

— = AF100475 Bayer et al., 
2000 


AF456795 — = Bayer & 
Cross, 2003 
EU385392 EU385107 EU385107 Panero & 
Funk, 2008 
EU749421 — - Fazekas 
et al., 2008 
= = U82019 Bayer & 
Starr, 1998 
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evolving genes of the chloroplast genome (Wolfe, 1991), is frequently employed for 
inferring plant phylogeny at the intergeneric level and above. Within Asteraceae 
matK, trnL intron, and trnL/F spacer have been used to infer higher level 
phylogenies, for example the trnL intron and trnL/F spacer at the tribal level (Bayer 
& Starr, 1998), matK for South African Gnaphalieae (Bayer et al., 2000), and matK, 
trnL intron, and trnL/F spacer for Australian Gnaphalieae (Bayer et al., 2002). We 
employed these plastid markers to reveal intergeneric relationships involving 
Leontopodium and relatives. 

The aims of this molecular study of Leontopodium, therefore, were to (1) elucidate 
phylogenetic relationships among species within the genus (including monophyly 
and infrageneric classification), (2) clarify relationships with other genera, including 
the recently described monotypic Sinoleontopodium Y.L.Chen (Chen, 1985; Chen & 
Yang, 2009), and (3) address the origin and evolutionary history of the alpine 
edelweiss (L. alpinum), a prominent taxon of the European mountain flora. 


MATERIALS AND METHODS 
Sampling 


We have included in our study 22 species of Leontopodium in 65 accessions (Table 1) 
from a broad geographic range and with a good representation from the Himalayan/ 
Tibetan centre of diversity. Less material was available from North and East Asia. 
Since many taxa occur only in remote areas, most DNA samples were obtained from 
herbarium specimens. Silica gel dried material from natural populations was 
available only from the European taxa (Leontopodium alpinum, L. nivale) and from 
the Taiwanese L. microphyllum. Vouchers are deposited in B and WU. Additional 
sequences, mostly of outgroup taxa, were drawn from the NCBI GenBank (Table 1). 

Several outgroup taxa were included to investigate the position of Leontopodium 
within Gnaphalieae: Anaphalioides subrigida, Anaphalis javanica, A. margaritacea, A. 
nubigena, Antennaria anaphaloides, A. carpatica, A. dioica, A. luzuloides, Castroviejoa 
frigida, C. montelinasana, Ewartia argentifolia, E. catipes, E. meredithiae, E 
planchonii, Ewartiothamnus sinclairii, Filago pyramidata, Gamochaeta pensylvanica, 
G. spicata, Gnaphalium hoppeanum, G. norvegicum, G. supinum, G. sylvaticum, 
Helichrysum leucopsideum, Leucogenes grandiceps, L. leontopodium and Pseudogna- 
phalium luteoalbum. In addition, several accessions from the ‘Relhania clade’ 
(Gnaphalieae) and from other Asteraceae tribes were added to root the plastid 
phylogenies, namely Leysera gnaphalodes, Oedera squarrosa and Relhania calycina 
(Relhania clade, Gnaphalieae), Achillea millefolium (Anthemideae), Solidago rugosa 
and Symphyotrichum novae-angliae (Astereae). 


DNA extraction and amplification 


Total genomic DNA was extracted from one leaf of herbarium material with the 
modified CTAB protocol of Doyle & Dickson (1987). The complete ITS (ITS1, 5.88 
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rRNA gene, ITS2), the 3’ part of ETS, a partial sequence of the plastid gene matK 
and the trnL/F region (trnL 5'-exon, intron, 3'-exon, intergenic spacer and trnF 
gene) were amplified with the universal primers ITS4, ITS5, ITS2, ITS3 (White et al., 
1990), 17SE and 26SE (Sun er al., 1994) for ITS; ETS-Hel-1 and 18SR (Baldwin & 
Markos, 1998) for ETS; trnL-c and trnL-c-1F/trnL-F (Taberlet et al., 1991) for trnL 
(UAA) intron and intergenic spacer between trnL (UAA), as well as 3'-exon and 
trnF (GAA) and 800F/1710R (Samuel er al., 2005) for a partial matK. Polymerase 
chain reaction was carried out in a volume of 25 ul, containing 10-50 ng of genomic 
DNA, 10 pmol of each primer, 1X Reddy Mix TM PCR Master Mix (including 
2.5 mM MgCl,; Abgene, Vienna) and 0.04% dimethyl sulfoxide (DMSO; for ITS and 
ETS) or 0.016% bovine serum albumin (BSA; for trnL/F and matK). The PCR 
profile consisted of an initial pre-melt at 94°C and 35 cycles of 1 min denaturing at 
94°C, 30 s annealing at 48°C and 1 min extension at 72°C followed by a final 
extension of 10 min at 72°C. Amplified fragments were checked on 1% agarose gel. 


Sequencing 


The amplification product was gel purified using Invitek Invisorb& Spin DNA 
Extraction Kit (Invitek Biotechnology & Biodesign Ltd, Berlin) The purified 
fragments were directly sequenced using dye terminator chemistry following the 
manufacturer's protocol and run on an ABI 377 automated sequencer (Applied 
Biosystems, Foster City, CA). The cycle sequencing reactions were performed with 
the same primers as used for PCR amplification. The programs Sequence Navigator 
version 1.0.1 and AutoAssembler version 1.4.0 (Applied Biosystems, Foster City, 
CA) were used to edit and assemble the complementary sequences. Alignments were 
obtained using the program Clustal X 1.5b (Thompson et al., 1997) and improved by 
manual refinement using the program BioEdit 7.09 (Hall, 1999). Pairwise genetic 
distances were calculated in the program BioEdit (Hall, 1999) using the feature 
Identity Matrix. Sequences are deposited in GenBank (www.ncbi.nlm.nih.gov). 


Phylogenetic analyses 


Fitch parsimony analyses (Fitch, 1971) were performed using PAUP* 4.0b10 
(Swofford, 2002). Four datasets were analysed: (1) ITS, (2) ETS, (3) partial sequence 
of mark, and (4) the trnL intron sequence and the trnL/F intergenic spacer. Heuristic 
search was conducted with equal weights, 1000 replicates of random sequence 
addition, tree bisection reconnection (TBR) branch swapping, and MulTrees on 
permitting 10 trees to be held in each step. Indels were treated as missing data, except 
for the ETS dataset where an additional parsimony analysis was conducted adding 
characters from gap coding using the MCIC method (Müller, 2006) as implemented 
in SequState (Müller, 2005). Trees were rooted in the plastid datasets using 
outgroups from the Gnaphalieae * Relhania clade’ (Leysera L., Oedera L., Relhania 
L'Hér.; Ward et al., 2009) and from more distantly related tribes, namely Achillea L. 
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(Anthemideae) and Solidago L. (in the matK dataset) or Symphyotrichum Nees (in 
the trnL/F dataset; both Astereae). The root for the nuclear markers was chosen 
using results of the plastid trnL/F dataset, in which Ewartia Beauverd was not sister 
to Leontopodium. To test the monophyly of Leontopodium, the following putatively 
close relatives from Gnaphalieae (according to Merxmiiller et al., 1977; Anderberg, 
1991; Bayer et al., 1996; Glenny & Wagstaff, 1997; Breitwieser et al., 1999) were 
included: Anaphalis DC. (ITS, ETS), Anaphalioides (ITS), Antennaria (ITS, ETS, 
matK, trnL/F), Castroviejoa Galbany, L.Sáez & Benedi (a recently recognised 
segregate of Helichrysum; ITS), Ewartia (ITS, ETS, matK, trnL/F), Ewartiothamnus 
Anderb. (ITS), Filago L. (ITS), Gamochaeta Wedd. (ITS, trnL/F), Gnaphalium L. 
(ITS, ETS), Helichrysum (ETS), Leucogenes Beauverd (ITS, trnL/F), and Pseudo- 
gnaphalium Kirp. (ITS). Nodal support was assessed via bootstrap values (Felsenstein, 
1985), which were calculated using PAUP* 4.0b10 with 10,000 bootstrap replicates 
each with 20 random sequence addition replicates holding a maximum of 10 trees per 
replicate, SPR branch swapping, and MulTrees on. 

To find the best-fit substitution models for the Bayesian analysis, all four datasets 
were tested using MrModeltest 2.2 (Nylander, 2005). The following partitionings 
were used for the analysis: (1) the ITS region was split into two partitions (the genic 
and the spacer regions of the ribosomal cistron) with GTR+I and GTR+G 
substitution models; (2) the matK/trnK region was also split into two partitions 
(the genic partial mark and the partial intron of the trnK gene) with GTR+G 
and F81 substitution models; (3) the trnL/F region was split into three partitions 
(the trnL intron, the trnL gene and the trnL/F spacer) with F81+G, the Jukes- 
Cantor, and GTR+G substitution models. For the ETS dataset Bayesian analysis 
was performed using the HKY+G with Mr Bayes (Ronquist & Huelsenbeck, 
2003), running four simultaneous Markov chains for 10,000,000 generations using 
four chains. Trees were sampled after every 1000 generations with a total of 40,000 
trees. The first 10% of trees was discarded to ensure that sufficient allowance was 
made for the chains to become stationary, thereby reducing the number of trees to 
36,000. 


RESULTS 
Nuclear markers 


The data matrix of the combined ITS region (ITS1, 5.8S, ITS2) included 22 ingroup 
taxa in 53 accessions (including one accession of Sinoleontopodium) and 23 outgroup 
accessions from 11 currently accepted genera (see Table 2 for analysis statistics and 
Table 3 for sequence statistics). Pairwise genetic distance in the complete ITS region 
(ITS1, 5.85, ITS2) ranged from 0.2% to 6.8% among the Leontopodium species. 
Pairwise sequence divergences between ingroup and outgroup were between 6.2% 
and 17%. Maximum parsimony and Bayesian analysis did show congruent results, 
with the latter providing considerably more resolved nodes. Therefore, phylogenetic 
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TABLE 2. Sequence statistics for ITS, ETS, matK and trnL/F 


Var. char. PICs No. of 
Marker Length (%) (%) CI RI trees Tree length 
ITS 649 30.8 21.4 0.68 0.88 710 381 
ETS 400/415? 35.8/38.3*  13/13.5*  0.9/0.91* 0.86/0.87* 9540/9200" 181/202? 
Partial 829-832 1.8/8.1" 0.84/3.6^ 1/0.95 Un. Oh 11/31° 15/795 
matK 
trnL 792 2.6/13.4^  0.5/3.5^ 10.94  1/0.91^ 1/85 22/125^ 
intron, 
trnL/F 


spacer 


“Including coded gap characters. 

"Including outgroups from tribes Anthemideae, Astereae and Gnaphalieae (Relhania clade). 
Abbreviations: var. char. = number of variable characters; PICs = number of parsimonious informative 
characters; CI = Consistency Index; RI = Retention Index. 


results of the ITS region are presented showing the Bayesian majority rule consensus 
tree (Fig. 1). 

Leontopodium was supported as monophyletic, with the inclusion of Sinoleonto- 
podium lingianum, the only species in Sinoleontopodium, with 100% bootstrap 
support (BS) and 1.00 posterior probability (PP) (Fig. 1). Some weakly 
(PP = 0.70) supported clades could be identified (clades A-C, Fig. 1), although 
the ingroup was not well resolved. Sister to the ingroup (BS 60, PP 0.91) were 
Antennaria, Castroviejoa, Gnaphalium, Filago and Gamochaeta (OG 1). Sister (BS 87, 
PP 0.98) to Leontopodium and the taxa of OG I were Anaphalis and Pseudognapha- 
lium luteoalbum, which together formed OG II (BS 100/PP 1.00). A third group of 
outgroup taxa comprised the Australian and New Zealand Gnaphalieae: Anapha- 
lioides, Ewartia, Ewartiothamnus and Leucogenes (BS 87, PP 0.98; OG II). 

The dataset of the ETS consisted of 39 ingroup accessions of 22 species and six 
outgroup accessions of five genera, namely Anaphalis, Antennaria, Gnaphalium, 
Ewartia and Helichrysum (Fig. 2). Generally, sequence divergence within the ETS 
dataset (Tables 2, 3) was extremely low except for Helichrysum and Ewartia. 
Antennaria dioica and Gnaphalium sylvaticum were nested within the ingroup. One 
clade (BS 57, PP 0.91) shared almost identical sequences and included Antennaria 


TABLE 3. ITS and ETS sequence characteristics showing length of the different ITS and ETS 
regions and GC content for ingroup and outgroup 


Length Length Length GC GC GC 


Region (aligned) ingroup outgroup mean (%) ingroup (%) outgroup (%) 
ITS1 257 252 251-254 43.7 42.7-45.5 41.5-45.5 
5.85 254 254 254 52.2 52.5 52.16 

ITS2 228 214217 214-223 50.9 50.23-51.63 45.2-53.6 


ETS 408 260-390 377-389 43.2 42-43.8 42.7-43.9 
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Fic. 1. Phylogenetic relationships of species of Leontopodium and related genera of tribe 
Gnaphalieae inferred from Bayesian and maximum parsimony analysis of the nuclear ITS 
region. Branches collapsing in the strict consensus tree of 710 equally most parsimonious trees 
are indicated by arrowheads. Bootstrap values/posterior probabilities are indicated at the 
nodes. Clades discussed in the text are marked with A-C and OG I-OG III. 


dioica, Gnaphalium sylvaticum, as well as most Leontopodium species (L. alpinum, L. 
artemisüfolium, L. franchetii, L. japonicum, L. leontopodinum, L. leontopodioides, L. 
microphyllum, L. nivale, L. ochroleucum, L. pusillum, L. sinense) and Sinoleontopo- 
dium lingianum. Two groups that were supported in the ITS analysis could be 
identified in the ETS analysis as well: Group B' (BS 74, PP 1.00) showed almost 
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Fic. 2. Phylogenetic relationships of species of Leontopodium and related genera of tribe 
Gnaphalieae inferred from Bayesian and maximum parsimony analysis of the nuclear 3’ ETS 
region. Branches collapsing in the strict consensus tree of 9200 equally most parsimonious 
trees are indicated by arrowheads. Bootstrap values/posterior probabilities are indicated at the 
nodes. Clades discussed in the text are marked with B’ and C. 
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identical sequences and included Leontopodium calocephalum, L. ochroleucum, L. 
souliei and L. jacotianum (the last came out in a different position in the ITS 
analysis). Group C contained Leontopodium sinense and L. dedekensii (BS 75, PP 
1.00) which share almost identical sequences. Leontopodium caespitosum appeared as 
sister to the ingroup comprised of Leontopodium, Antennaria and Gnaphalium (which 
were supported with BS 60, PP 0.92). The ingroup was supported with BS 100, PP 
1.00, when rooted with Ewartia, having Helichrysum leucopsideum, Ewartia catipes 
and Anaphalis nubigena as clearly defined outgroups. 


Plastid markers 


Both plastid markers used in this study (matK, trnL intron—trnL/F spacer) had very 
low sequence divergences, with ingroup sequences being identical or almost identical. 
The matK dataset comprised 17 accessions of the ingroup and two outgroup 
accessions from the Gnaphalieae, namely Antennaria luzuloides and Ewartia catipes 
(for analysis details refer to Table 2). Furthermore, five outgroup accessions from 
other tribes of Asteraceae were included. 

The two outgroup genera from Gnaphalieae were clearly separated by three 
synapomorphic substitutions as compared to Leontopodium. Considering the low 
variation within the whole dataset, these substitutions provide further support for 
monophyly of the ingroup (BS 77, PP 0.98; Fig. 3). 

Within Leontopodium, one group consisting of L. cf. stracheyi (3), L. ochroleucum 
(1 and 2), L. andersonii (4), and L. leontopodioides (3) was supported (BS 64, PP 
0.98). However, the relevance of this group might be in doubt as they share only 
a single nucleotide substitution not present in other Leontopodium accessions. This 
well-supported group, which is not supported by either morphology or chorology, 
may be an artefact of the low sequence divergence. Leontopodium microphyllum 
appears sister to the remainder of Leontopodium (BS 52, PP 0.88). 

The combined dataset of the trnL intron and the trnL/F spacer region included 10 
ingroup accessions and three additional accessions from the crown group of the 
Gnaphalieae (Ward et al., 2009), namely Antennaria carpatica, Ewartia catipes and 
Leucogenes grandiceps. The trees were rooted using three accessions from the 
* Relhania clade’ of Gnaphalieae (Ward et al., 2009) and one accession each from 
Astereae and Anthemideae (analysis details in Table 2). All accessions of the ingroup 
had identical sequences for the intron, exon and spacer regions, except for 
Leontopodium microphyllum, which showed two autapomorphic substitutions and 
one substitution shared with the other included genera of Gnaphalieae. The ingroup 
was supported by BS 67 and PP 0.99. Sister to Leontopodium was Antennaria dioica 
(BS 68, PP 0.97; Fig. 4). Ewartia, Gamochaeta and Leucogenes were together sister to 
Antennaria and Leontopodium (BS 97, PP 1.00). The outgroup taxa of Gnaphalieae 
shared one synapomorphic substitution. 

Despite the lack of resolution, the plastid sequences support the monophyly of 
Leontopodium. Furthermore, the sequences from Leontopodium microphyllum were 


MOLECULAR PHYLOGENY OF THE EDELWEISS 253 





. cf. stracheyi 3 
57/0.99 . ochroleucum 1 
64/0.98 . andersonii 4 

. leontopodioides 3 
. ochroleucum 2 

. alpinum 11 

. caespitosum 2 

. calocephalum 2 
52/0.78 nanum 2 


nivale 1 


WNIAOdOLNO37 


. ochroleucum 3 
. pusillum 3 
. sinense 2 


3V3lTVHdvN9 


77/0.98 


souliei 
. stracheyi 1 
97/1.00 . microphyllum 1 


. microphyllum 2 


52/0.88 
« 


64/1.00 Antennaria luzuloides 


rr FE FEEF FEEF FEEF E 


Ewartia catipes 1 
>50/0.61 Relhania calycina 
78/.00 «+ Leysera gnaphalodes 2 
Oedera squarrosa 2 
Achillea millefolium 2 E 
Solidago rugosa H 


Fic. 3. Phylogenetic relationships of species of Leontopodium and related genera of tribe 
Gnaphalieae, as well as of additional accessions from the tribes Astereae and Anthemideae 
inferred from Bayesian and maximum parsimony analysis of the partial plastid matK region. 
Branches collapsing in the strict consensus tree of 31 equally most parsimonious trees are 
indicated by arrowheads. Bootstrap values/posterior probabilities are indicated at the nodes. 
AST = Astereae, ANT = Anthemideae. 


distinct from the other sequences of the ingroup in both regions. Further infrageneric 
grouping could not be inferred using these markers. However, it is noteworthy that 
the ingroup sequences were clearly different from the outgroup sequences of 
Antennaria, Ewartia, Gamochaeta and Leucogenes and that according to the trnL/F 
dataset, Antennaria was the closest relative of Leontopodium of the four included 
genera of the crown radiation of Gnaphalieae (according to Ward et al., 2009). 


DISCUSSION 
Intergeneric relationships 


The genus Leontopodium is monophyletic if the monospecific Sinoleontopodium is 
included. Three of the four employed markers support its monophyly: ITS (BS 100, 
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Fic. 4. Phylogenetic relationships of species of Leontopodium and related genera of tribe 
Gnaphalieae, as well as of additional accessions from the tribes Astereae and Anthemideae 
inferred from Bayesian and maximum parsimony analysis of the plastid trnL gene and trnL/F 
spacer region. Branches collapsing in the strict consensus tree of 8 equally most parsimonious 
trees are indicated by arrowheads. Bootstrap values/posterior probabilities are indicated at the 
nodes. AST — Astereae, ANT — Anthemideae. 


PP 1.00), matK (BS 77, PP 0.98) and trnL/F (BS 67, PP 0.99). The ETS dataset 
included part of Antennaria and Gnaphalium in the ingroup (BS 100, PP 1.00), 
probably due to the lack of sequence divergence and homoplasy. 

The position of Sinoleontopodium within Leontopodium is morphologically and 
chorologically supported and, therefore, the following new combination is made: 


Leontopodium lingianum (Y.L.Chen) Dickoré, comb. nov. — Sinoleontopodium 
lingianum Y .L.Chen, Novon 19: 24 (2009) [Sinoleontopodium lingianum Y .L.Chen, 
Acta Phytotaxon. Sin. 23: 458, fig. 1 (1985), nom. inval. (two types cited)]. 


Leontopodium lingianum is a local endemic of the high mountains around the 
Tsangpo Gorge of Southeast Tibet (Mainling, Mt Namcha Barwa, Mt Gyala Peri), 
where it occurs at 4400-4850 m in a perhumid alpine belt. Leontopodium lingianum 
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displays distinct morphological adaptations: a dense cushion of columnar shoots; 
leaves many, imbricate, small, ovate-spathulate, long and densely lanate; calathidia 
sessile, solitary, without a star. However, all individual characters and its distribu- 
tion are contiguous to other species of Leontopodium. The peculiar habit of 
Leontopodium lingianum is closely matched by L. aurantiacum, a species from the 
Burma/Yunnan border, which was not available for the present study. Complete 
dioecy, one of the main distinguishing characters of Sinoleontopodium, occurs 
facultatively in Leontopodium haastioides, another acaulescent, cushion-forming 
species from the Central Himalayas, and in several other species of Leontopodium. 
Hirsute corollas, although not to the extent of Leontopodium lingianum, are also 
found in other species of Leontopodium. The exact position of Leontopodium 
lingianum within the genus is still unclear; ITS1 data place it within a large polytomy. 

Relationships of Leontopodium within Gnaphalieae are currently largely un- 
resolved. Based on cladistic analysis of morphological data, Anderberg (1991) 
concluded that the ‘New Zealand edelweiss’ Leucogenes is sister to Leontopodium. 
This hypothesis was, however, rejected by the molecular phylogenetic analyses of 
Glenny & Wagstaff (1997) and Breitwieser et al. (1999). Based on morphology 
Merxmiiller et al. (1977) placed Leontopodium within an * Anaphalis-group' consisting 
of Antennaria, Anaphalis and Leontopodium. However, it grouped with Antennaria in 
the molecular study of Glenny & Wagstaff (1997) as sister to Pseudognaphalium and 
Anaphalis. The latter findings are partly concordant with Bayer et al. (1996) who 
identified Antennaria as sister to Leontopodium and both as sister to the Australian— 
New Zealand genus Ewartia. The studies of Ward er al. (2009) on the phylogeny of 
Gnaphalieae, with a focus on African and Australasian taxa, show that relation- 
ships within Gnaphalieae are not well resolved, that taxa from the so-called Relhania 
clade are basal to the rest of the Gnaphalieae, and that Leontopodium is sister to 
Antennaria. Our combined plastid and nuclear data indicate that Leontopodium is 
closely related to a group comprised of representatives from the bicentric (south and 
north hemisphere) genus Gnaphalium, the northern hemisphere genus Antennaria, 
the newly described genus Castroviejoa (Galbany-Casals et al., 2004b), and Filago. 
Accordingly, none of the solely southern hemisphere (mainly Australian-New 
Zealand) taxa, which form a distinct clade in our study and include Ewartia which 
was sister to Leontopodium and Antennaria in the study of Bayer et al. (1996), is sister 
to Leontopodium. Despite strong morphological similarities with Leontopodium, the 
large genus Anaphalis, primarily of temperate East Asian—North American distri- 
bution but also with a considerable representation in tropical Indo-Malaysia, seems 
to be even more distantly related to Leontopodium and instead is rather close to 
Pseudognaphalium. 


Species delimitation and infrageneric classification 


The ITS dataset, and to a lesser degree the other markers, reveals information about 
infrageneric relationships within Leontopodium. Three groups may be recognised but 
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these do not correspond well to the hierarchical scheme of Handel-Mazzetti (1927). 
While species delimitations within Leontopodium were partly supported by our study 
(not all taxa were monophyletic), broader infrageneric relationships were not fully 
resolved, indicating close species relationships and suggesting recent speciation. 
Hybridisation seems probable in Leontopodium, due to nested distributions of 
species and often contiguous variation, especially in the Himalayas and the 
Hengduan Shan. We could not confirm any of the hybrids proposed by Handel- 
Mazzetti (1927). Furthermore, species delimitation in Leontopodium seems to be 
complicated by the possible occurrence of apomixis. This was mentioned, for 
example, in the Flora of Pakistan (Qaiser & Abid, 2003). The only species where 
apomixis has been studied and found is in Leontopodium alpinum (Sokolowska- 
Kulczycka, 1959; Maugini, 1962). In a review of apomixis in Asteraceae, Noyes 
(2007) mentions that apomixis within Leontopodium has been poorly studied so far. 
On the basis of available data, therefore, no further conclusion on apomixis can be 
drawn. 

We have identified the following molecular clades of Leontopodium, each with 
morphological integrity: 


Group A contains the European taxa (Leontopodium alpinum, L. nivale). These taxa 
belong to the type section Leontopodium [Alpina Hand.-Mazz.]. They are genetically 
distinct from all Asian species in the ITS analysis, although perhaps not as 
distinctly separated as might be expected from the geographic distance. In ETS, 
matK and combined trnL intron-trnL/F spacer datasets, sequences of the European 
Leontopodium alpinum were the same as for many Asian Leontopodium species. In 
contrast, the similarly geographically isolated Leontopodium microphyllum from 
Taiwan had clearly distinct sequences in all three datasets. The question is still 
unresolved as to whether the European taxa comprise two distinct species (Handel- 
Mazzetti, 1927), subspecies (Tutin, 1973; Greuter, 2003), or a series of varieties or 
unclassifiable forms. We recommend recognising Leontopodium alpinum and L. 
nivale at species rank until more sophisticated population-level analyses may clarify 
their status (S. Safer et al., unpubl.). Similar problems also occur in Asian 
Leontopodium. 


Group B, consisting of Leontopodium calocephalum, L. leontopodioides, L. ochroleu- 
cum and L. souliei (all in Leontopodium sect. Leontopodium in Handel-Mazzetti, 
1927), is moderately supported by ITS and ETS, in the latter with the inclusion of L. 
jacotianum. Morphologically, this group (except Leontopodium jacotianum) is quite 
coherent and close to the European representatives. Except for Leontopodium 
jacotianum, group B species occur mostly in areas with cold, humid or mesic 
conditions in Central and Northern Asia and are largely absent from the warm and 
humid to perhumid areas of the Sino-Himalayan region and the dry-cold core of 
Central Asia. Morphological, ecological and chorological properties would thus 
make group B species (except for Leontopodium jacotianum) candidates for the 
closest Asian relatives of the European taxa. 
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Group C, supported by both the ITS and ETS datasets, consists of Leontopodium 
sinense and L. dedekensii. Both species belong to Leontopodium sect. Nobilia 
(Handel-Mazzetti, 1927) and are very similar in habit and range, being partly 
sympatric, from subtropical Southwest China to subalpine altitudes of the Eastern 
Tibetan Plateau. Although Leontopodium sinense and L. dedekensii are morpholog- 
ically and ecologically sufficiently distinct, they share identical sequences in all the 
markers used in this study. Although Handel-Mazzetti’s Leontopodium sect. Nobilia 
is morphologically well defined, the remaining species of the section are in 
unsupported positions in both the ETS and ITS phylogeny. 


Origin and evolution 


In the absence of any fossil record, a hypothesis for the origin of Leontopodium needs 
to take into account our molecular findings plus inferences from present-day 
ecology, distribution and diversification. In the molecular data two morphologically 
similar species, Leontopodium caespitosum (from the Hengduan Shan) and L. 
microphyllum (Taiwan), were basal, suggesting they are part of an ancestral stock 
within the genus. The East Himalaya-Hengduan Shan region on the south and east 
rim of the Tibetan Plateau, home of Leontopodium caespitosum, comprises a major 
hotspot of biodiversity outside the tropics and is renowned for the occurrence of 
numerous ‘relict’ taxa (Liu et al, 2006). The island flora of Taiwan, where 
Leontopodium microphyllum occurs, is dominated by paleotropical elements, while 
its mountains show strong biogeographic parallels to the Himalaya-Hengduan Shan 
region (Wu & Wu, 1998). The southern edge of the temperate Sino-Japanese 
(Grisebach, 1872; Good, 1974) or East Asiatic floristic region (Diels, 1901; 
Takhtajan, 1969) towards the paleotropics of Indo-Malesia includes the highest 
present-day species diversity in Leontopodium and may also have been its centre of 
origin. 

Many species of Leontopodium show distinct ecological adaptations to high 
altitudes and dry-cold conditions, but this may not have been the climate of its 
origin. In its Southwest Chinese centre of diversity, the genus spans a wide range of 
elevations and climates. The ecology of Leontopodium shows strong links to what is 
commonly called an ‘Arcto-Tertiary laurel forest flora’ (Takhtajan, 1969), and the 
highest concentration of species occurs at or near the subtropical/meridional laurel 
forest region of Southwest China. The unique floristic richness of this region can be 
explained by an unbroken connection between tropical and temperate, generally 
rather favourable, climates, with wide altitudinal gradients and relative stability 
throughout a long geological timespan (Axelrod et al., 1998; Kubitzki & Krutzsch, 
1998). The modern stock of Leontopodium may have evolved from plants adapted to 
an equable, subtropical climate, which prevails in part of this area. 

Speciation of Leontopodium has been influenced by the turbulent geological 
history in its centre of distribution, the Tibetan Plateau. Following the collision of 
the Indian Plate with Eurasia, the Tibetan Plateau began to uplift between 70 Ma 
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(Unsworth et al., 2005) and 50 Ma ago (Searle, 1991; Royden et al., 2008). The uplift 
of the Himalayas and the Tibetan Plateau led to large-scale climate changes in Asia, 
basically by establishing the monsoon system, but also causing a general transition 
towards cooler, more pronounced seasonal climates over large areas (Jen, 1982; 
Axelrod et al., 1998; Kubitzki & Krutzsch, 1998). During the Pleistocene several 
glaciations affected the region (Zhou et al., 2006). While the (‘Sino-Japanese’) flora 
of the Hengduan Shan is considered to be derived from the ancient Laurasian flora, 
epibiontic and relict in nature, the Himalayan (‘Sino-Himalayan’) flora may be 
characterised as being mixed with neoendemic elements derived from several 
(Laurasian, Tethyan and Gondwanic) geoelements (Wu & Wu, 1998). Species 
diversity in Leontopodium is higher on the Hengduan Shan eastern edge of the 
Tibetan Plateau, and individual species’ distribution patterns would also imply 
a subsequent colonisation of the Himalayas, approximately separated by the 
Mekong-Salween divide (Ward, 1921). However, a considerable distribution gap 
between the regions, such as in the genus Roscoea Sm. (Ngamriabsakul et al., 2000), 
is absent in Leontopodium, possibly because this genus spread much more into the 
cold-arid hinterland. Distribution and diversification of Leontopodium aligns to 
climatic and geomorphological patterns around the Tibetan Plateau, and repeated, 
partial re-colonisations of large high altitude areas (ie. during the Pleistocene 
glaciation cycles) would seem most likely. 

Recent molecular phylogenetic studies in various Asian alpine plant groups have 
also assumed that extreme abiotic factors related to the geological and climatic 
history of this region have substantially affected biota and speciation processes (Liu 
et al., 2002, 2006; Wang & Liu, 2004; Y. J. Wang et al., 2004, 2007; A. L. Wang et al., 
2005). In all these studies levels of nucleotide substitution and support for internal 
branches in the phylogenetic analyses are low, and morphologically diverse species 
often group together. These patterns could also be detected in this study. They have 
perhaps arisen due to rapid radiation and hybridisation triggered by multiple, severe 
climatic changes, and habitat fragmentation and rejoining. Ecological constraints 
might also have been responsible for the occurrence of similar morphological 
characters in unrelated lineages, resulting in the grouping of taxa with few 
morphological similarities. 

The origin of the European taxa of Leontopodium cannot be inferred with high 
confidence on the basis of available data. Nevertheless, a migration from Southeast 
Asia via Middle Asia seems most probable. While there are no extant intermediate 
occurrences between Central Asia and Eastern Europe, long-distance dispersal 
spanning this gap does not seem likely as the pappus of achenes in Leontopodium 
cypselas is not well suited for wind dispersal. The geological history of Eurasia and 
the alpine orogenic belt, however, might suggest migration during cooler climatic 
cycles of the Pleistocene. A continuous distribution from Central Asia to Eastern 
Europe might have been possible during the last interglacial (120,000 years ago) 
along with widespread herb-grass and Artemisia-grass steppe formations (Grichuk, 
1992). 
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